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We investigate the Josephson effects in the junction formed by Z3///-class topological and s-wave 
superconductors, by embedding a quantum dot in the junction. Three dot-superconductor coupling 
manners are considered, respectively. As a result, the Josephson current is found to oscillate in 27r 
period. Moreover, the presence of Majorana doublet in the Dlll-class superconductor renders the 
current finite at the case of zero phase difference, with its sign determined by the fermion parity of 
such a jucntion. In addition, the dot-superconductor coupling plays a nontrivial role in adjusting 
the Josephson current. When the s-wave superconductor couples to the dot in the weak limit, the 
current direction will have an opportunity to reverse. The results in this work will be helpful for 
understanding the transport properties of the Dill-class superconductor. 
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I. INTRODUCTION 

Topological superconductors (TSs) have received much attentions from experimental and theoretical as¬ 
pects because Majorana modes appear at the ends of the one-dimensional TS which can potentially be used 
for topological quantum computation.^ilf^ Due to the possibility of achieving Majorana modes, the systems 
with TSs show abundant and interesting physical characteristics .1^ For instance, in the proximity-coupled 
semiconductor-TS devices, the Majorana zero modes induce the zero-bias anomaly.^ A more compelling TS 
signature is the unusual Josephson current-phase relation. Namely, when the normal s-wave superconductor 
nano-wire is replaced by a TS wire with Majorana zero modes, the current-phase relation will be modified to 
be Ij ~ V sin ^ with the period Jtt {(j) is the superconducting phase difference and V is the fermion parity). 
This is the so-called the fractional Josephson effect.l^^^ 

More recently, the time-reversal invariant TSs, i.e., the Dill symmetry-class TSs,l^®ti3 have become one 
new concern.lJSHIZl jn such a kind of TSs, the zero modes appear in pairs due to Kramers’s theorem, dif¬ 
ferently from the chiral TSs. Gonsequently, two Majorana bound states (MBSs) will b e localized at each 
end of the DIII-TS nanowire, leading to the formation of one Kramers doublet.l^^E^ Since the Kramers 
doublet is protected by the time-reversal symmetry, it will inevitably drive some new transport phenomena, 
opposite to the single Majorana zero mode. Up to now, many groups have proposed proposals to achieve 
the DIII-TS nanowires, by using the proximity effects of d-wave, p-wave, s±-wave, or conventional s-wave 
superconductors.lSSHlIl Meanwhile, physicists have paid attention to the quantum transport phenomena con¬ 
tributed by the Kramers doublet, and some important results have been reported.l^^^^ For instance, in the 
Josephson junction formed by the Majorana doublet, the period of the Josephson current will be varied by 
the change of fermion parity in this system.l^ This exactly embodies the nontrivial role of Majorana dou¬ 
blet in contributing to the Josephson effect. However, for further understanding the property of Majorana 
doublet, Josephson effects in any new junctions should be investigated. 

In this work, we would like to discuss the Josephson effect in a complicate junction, i.e., one junction 
formed by the indirect coupling between a Dill-class TS and s-wave superconductor with one embedded 
quantum dot (QD). Our calculations show that the interplay between the Dlll-class TS and the s-wave 
superconductor induce interesting results. Although the Josephson current oscillates in 27r period, the 
presence of Majorana doublet in the Dlll-class TS renders it finite at the case of zero phase difference, with 
its sign determined by the fermion parity of the whole system. Also, QD plays a nontrivial role in adjusting 
the Josephson effect. To be concrete, in the extreme case where the s-wave superconductor couples weakly 
to the QD, the direction of the Josephson current will be reversed. All these results embodies the specihc 
role of Dlll-class TS in driving the Josephson effect. 



FIG. 1: The junction formed by Dlll-cless topological and s-wave superconductors. One QD is embedded in the 
junction. 
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II. THEORETICAL MODEL 


The considered Josephson junction is illustrated in Figj^ where one _D///-class TS couples to one s-wave 
superconductor via one QD. The Hamiltonian of such a system can be written as H = Hp + Hs + Hjj + Ht- 
The first two terms, i.e., Hp and Hs, denote the Hamiltonians of the DHI-c\ass TS and s-wave superconductor 
respectively, which is written aJ^ 

Hp = ^ ^ ~ ^ ^ + d~ h.c.) -f- ^ ^ [(~'io'ig2)erfT'd~ h.C.], 

ja ja jaa' 

Hs = (1) 

jcr ja j 


and /J^ {cja- and fj„) are the electron creation (annihilation) operators in the Dili-class TS and s-wave 
superconductor, respectively. /Xp and fis are the chemical potentials of the superconductors, and Ap and Ag 
are the Copper-pair hopping terms. Hd describes the Hamiltonian of the QD. For a single-level QD, it takes 
the form as Ho = Hdo + H^e with 

Hdc=eo^{na-^), Hee = U{n^ ( 2 ) 


ria = d\da is the electron-number operator, in which d\ and da are the creation and annihilation operators 
in the QD. Eq is the QD level, and U denotes the intradot Coulomb repulsion. Ht, the last term of H, 
represents the couplings between the QD and the superconductors. It can be given by 

Ht =+h.c.)+h-c-): (3) 


where It and ts are the QD-superconductor coupling coefficients, respectively. 

It is well-known that the phase difference between superconductors drives finite current through one 
Josephson junction. With respect to such a junction, the current properties can be evaluated by the following 
formula 


_ 2e d{H) 
~ li de 


(4) 


In this equation 9 = 9t — Os Is the phase difference between the superconductors, and (• • •) is the thermal 
average. As a typical case, i.e., the zero-temperature limit, the Josephson current can be simplified as 


_ 2e dEcs 

h do 


(5) 


This formula shows that the calculation about the Josephson current is dependent on the ground-state (GS) 
level of this system. 

In the low-energy region, Dlll-dass TS only contributes Majorana doublets to the Josephson effect, so 
Hp describes the coupling between two Majorana doublets. For the extreme case of one infinitely-long TS, 
the coupling strength between the Majorana doublets decreases to zero.® Following this idea, we project 
Ht onto its zero-energy subspace. As a result, H can be rewritten as 

H = £o - ^) + U{n^ - - ^) + fLfka +As fit f-K 

a ka k 

+tT '^{E^^^da - e~^^^^dl)-foa - ^ Vksflada + h.c., (6) 


where Hs has been projected into in the Bloch space with = —/x' — 2t' cos k and Vks = J2k loa, 

the Majorana operator, which obeys the anti-commutation relationship of {7o<T,7ocr'} = “^Aaa'- Based on the 
renewed expression of H, we next try to diagonalize the Hamiltonian of such a Josephson junction. 


III. DIAGONALIZATION OF THE JUNCTION HAMILTONIAN 

The continuum state in the s-wave superconductor hinders the diagonalization of the system’s Hamiltonian. 
In order to present a comprehensive analysis, we would like to consider three cases, i.e., the cases of It ^ ts, 
It ~ ts, and tT ts, followed by the utilization of different approximation methods. The following are the 
detailed discussion processes. For convenience, they are named as Case I, Case H, and Case HI, respectively. 
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A. Diagonalization of H in Case I 


In Case I where tr ^ ts, the subsystem formed by the QD and Majorana doublet can be considered to be 
one system, whereas the s-wave superconductor can be viewed as a perturbation factor. We next simplify 
the system Hamiltonian by means of the perturbation theory. Ignoring the Coulomb interaction term in the 
QD, we can write out the action of the subsystem of QD and s-wave superconductor 




(7) 


= 


and = 


. With the action S', we can express the partition function as a path integral. 


where Hdo = SqO'z, 'Hs = ^kO'z + ‘^sO'x, and VkS = Vkso'z- As for the field operators, they are given by 
fkt 

fU- 

i.e., Z = in which the measure P'l'fc denotes all the possible integral paths. 

After integrating out the fermion field '1'^, with a Gaussian integral, the partition function will become 
a “generating functional” 
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Zl'i’d, = det[9^ -I- Hdo] exp / drdr' J^k r - T')vlg'i>d{T') 


( 8 ) 


G{k, t — t') = —{Tr'i'k{T) 0 'I'j.(T')) is defined in the s-wave superconductor. It obeys the Fourier expansion 


G{k,T) = ^J2n^ with G(A:, = [det(iw„ - Hs)] ^ 

to write out the effective expression of the action in the Fourier space directly, i.e.. 


-I- ^k AsC' 


i0 


AgC 'l^n 


This allows us 


Seff = -^TrlnG-i(iw„) - f ^dG(/c, 


k.ioJn 


r. .2 , .21 -h 


(9) 


With the help of the expression of SeS, it is not difficult for us to get the perturbative Hamiltonian of the 
s-wave superconductor on the QD 


H' ~ 

N 


N 


/tr_ 

^ d\.c2 


Cfc 


e^ + A? 


A, 

?^ + a4 




( 10 ) 


Via a unitary transformation, the system Hamiltonian can be expressed as the following form 

HeS = so'^iria - 4 +As{e’-^d\dl + + U{n^- - 4 + tr'^id^ - 4 ) 7 o < t , (H) 


As 


2 " " ski/' \ I 2 2 

(7 

Such a result indicates that the weak Andreev reflection between the s-wave 


with A, = 

superconductor and QD induces a weak s-wave pairing potential on the QD, which is exactly the so-called 
the proximity effect. 

For the sake of diagonalizing such a Hamiltonian, we need to introduce local Majorana operators rjis 
through da = + t?? 2 cr) and dj. = ~ iV 2 a)- And then, by defining Dirac fermionic operators 

a = + ^ = ;^(’ 72 t+ *^ 2 ;), and / = ^( 7 oT + * 7 o;), we can obtain a new expression of ideff, i-e.. 


iJeff = As sind(a^a — 6^5)-I-(eo-I-iAg cos0)a^6-I -(£0 — tAg cos0)6^c 
-U{na - ^){nb ~ V2tT{b'^f + h.c.). 


( 12 ) 


According to Eq.(12), the BogoliubovCde Gennes equation iFeff^ = E'^ can be built, and then the 
eigenvalues of Heff can be worked out. On the basis of {|000), |001), |010), |100), |110), |101), |011), |111)}, 
the matrix form of H can be obtained (|nan.;,n/) = \na)\nh)\n, where Ua = a^a, rik = b^b, and — f^f)- 
Note that in the TS-existed system, only the parity of the average particle occupation number is the good 
quantum number, thus the matrix form of HeS should be given according to fermion parity (FP). In the 
case of odd FP, 'I'o = ci|001) -I- C 2 IOIO) -I- C 3 IIOO) -I- C 4 IIII), and then 


= 


-U ^-^tT 0 0 

— — AgSind-l-^ eo — iAsCos9 0 

0 eo+tAgCosd As sin 0-1-^ 0 

0 0 0 


(13) 
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In the case of even FP, 'I'e = ci|000) + C 2 IOII) + C 3 IIOI) + C 4 IIIO), and 


Ks’ = 


0 00 

0 — AsSin0+^ eo — iAsCos6 0 

0 £o + *^sCos0 AsSin0+^ — 

0 0 -V2tT 


(14) 


It is easy to find that (9) = (0 + tt). Thus, the Josephson effect can be clarified by only analyzing 

the current oscillation result in one FP. 


B. Diagonalization of H in Case II 


In Case II where tx ^ ts, H is difficult to diagonalize due to the presence of continuum state in the s-wave 
superconductor. However, according to the previous works, the zero band-width approximation is feasible 
to solve the Josephson effect contributed by the s-wave superconductor.^^ Within such an approximation, 
the Hamiltonian can be simplified as 


Hes = eo - i) + t/(nt - - ^) + ^ 

(7 (7 

+tT'^{da - dlhoa - ^ tsfld^ + h.C.. 


(15) 


By defining + and /] = _ j^^cr) with a = + and /3 = ^(C 2 t+< 2 ;), 

we get the Hamiltonian in the spinless-fermion representation 

HeS = £o{a^b + b'^a)-U{na-^){nb-^) 


-l-Ag sin0(a^a — /I’^/3) -I- (^ -I- iAg cos 9)a^/3 -|- (^ — iAg cos 9)f3^a 
—V2tx{b^f + pb) — tg{b^a + a^b + a) P + pa). 


On the basis of odd FP, the matrix of can be expressed as 




0 K 
0 


0 0 
(O) Q 

in ^ 

0 


where 


and 


''-III — 



r ^ 

4 

£0 


0 


-ts 


0 



£0 

u 

4 


-ts 


0 


— p2tT 


) ^ 

0 

-ts 

A« 

sin 0 

u 

4 

^ -1- iAs cos 9 

0 



-ts 

0 

C- 

iAs 

ZOS0 

-As sind 

u 

4 

0 



0 

-V^tj 


0 


0 


u 

4 J 


El 

In 

0 

0 

0 

p2tT 

-ts 

0 

0 

0 

t*N 

E 2 

0 

-ts 

0 

0 

0 

0 

0 

0 

0 

0 

u 

4 

0 

0 

-ts 

0 

0 

ts 

0 

0 

-ts 

0 

u 

4 

0 

0 

£0 

0 

0 

0 

0 

0 

0 

0 

-El 

t*N 

0 

0 

£0 

0 

pitT 0 

-ts 

0 

tN 

-E 2 

0 

£0 

0 

ts 

-ts 

0 

0 

£0 

0 

0 

u 

4 

0 

0 - 

\/2t 

0 

0 

0 

0 

0 

£0 

0 

-E 2 

tN 

0 

0 

0 

ts 

0 

£0 

0 

0 

tN 

-El 

-ts 

0 

0 

0 

0 

0 

ts 

— p2tT 

0 

-ts 

u 

4 


(16) 


(17) 
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with El = As sin6* — E 2 
of can be given by 


= —As sin0 — and In = ^ + *As cos 6*. On the basis of even FP, the matrix 


Hi^ = 


u 

4 

0 

0 


0 
0 


0 

0 

'(e) 

'■// 


where 


and 


"'ll — 



r ^ 

4 

£0 


0 

-ts 


0 




£0 

u 

4 


-ts 

0 


— \/2tT 



II 

0 

-ts 


-As sin 0 - F 

^ - iAs 

COS0 

0 




-ts 

0 


^ + iAg cos 0 

As sin0 

u 

4 

0 




0 

-y/2t 

T 

0 

0 


u 

4 



E2 

t*N 

0 

0 

0 

-ts 

0 

0 


0 

In 

El 

0 

-ts 

0 0 

0 

0 

0 


0 

0 

0 

u 

4 

0 

0 -ts 

0 

0 

ts 


0 

0 

-ts 

0 

u 

4 

0 0 

£0 

0 

0 


0 

0 

0 

0 

0 

— E2 tN 

0 

0 

£0 


0 

y/2tT 

0 

-ts 

0 

t*N -El 

0 

£0 

0 

ts 

-ts 

0 

0 

eo 

0 0 

u 

4 

0 

0 - 

y/2t 

0 

0 

0 

0 

0 £0 

0 

-El 

tN 


0 

0 

0 

ts 

0 

£0 0 

0 

tjn 

— E2 


ts 

0 

0 

0 

0 

0 ts 

-V^tr 

0 

-ts 


u 

4 


Therefore, the different-FP matrix forms of Heff obeys the relationship (0) = (0 - 

result in Case L 


r(e)/' 


(18) 


- tt), similar to the 


C. Diagonalization of H in Case III 


In Case III, we turn to the discussion about the diagonalization of ff when tx ^ts- In such a case, the 
QD will dip in the s-wave superconductor, leading to the formation of a composite s-wave superconductor. 
Consequently, the considered structure will be transformed into a junction in which the Majorana doublet 
couples to a s-wave superconductor directly. Its Hamiltonian can thus be written as 


- 4J7o. + E + As 


(19) 


k(7 


In Eq.(19l, F^a originates from the unitary transformation that da- = "kFka and fua = J2k' dkk'Fk'a, and 
Wk = Vktr denotes the coupling between the Majorana doublet and the s-wave superconductor. It is easy to 
prove that in composite s-wave superconductor couples weakly to the Majorana doublet [See the appendix]. 
As a result, the composite s-wave superconductor can be considered as perturbation. With respect to the 
Hamiltonian in Eq.(|I9[), the action can be written as 




in which TLs = 

a Ase*« ■ 
.Ase-*« -a _ 

and >V(fc) = VFfcCr^ with = 

1 1 

1 _1 

and Fq = 

1 1 

0 0 

derivation process in Case I, we 

express the partition function as 

a path 

integral 



( 20 ) 

. Similar to the 


Z = / = I V^kV'H' exp 


- / rfr E{^l[5r+Hs(fc)]^-fe + (WFj§fc+Wt§tro)} 

^ k 

Integrating out the fermion field with a Gaussian integral, we simplify the partition function as 

Z[ro,rJ,] = det[9^]exp j dTdT'^Tl[T)WG{k,T- t')W^Tq[t') 

k 

= det[clT-]exp J dTdT'rl{T)iTGd{T — T')tlpTo{T') 


( 21 ) 
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where G{k,T- t') = -{Tr^uir) ® ^fc(T')), Gdir - t') = -(T^^'d(r) ® and ir = Itu^. 

Since Gd('r) obeys the relationship that Gdir) = ^ J2n in the Fourier space the effective 

action can be expressed as S'eg = - Tr ln[a;^] - ^ rJ)(ia;„)fTGd(fc, Accordingly, 

the Josephson Hamiltonian in Case III can be given by Heg{6) = —^ 'l2iuj [2*IniGd,/ie(fa;„)] 7ot7o^- At the 
zero-frequency limit, the approximated form of H^g can be written as H^g = J(ijot^oi) sin0 with 

J = « As4/4- (22) 

By defining / = ^(7ot + *7o^)i we obtain the result that Hgg = J{nj— sin0, i.e., = (+/—)7 sin0. 

Therefore, the Josephson current can be directly written as ) J cos0. Surely, such a result 

is consistent with that in RefP^. 



FIG. 2: Spectra of odd-FP Josephson current in Case I of [7 = 0 and U = 2.0, respectively. The coupling strength 
between the QD and superconductors are taken to be It = 0.5 and = 0.04. (a)-(b) [7 = 0; (c)-(d) U = 2.0. 


IV. NUMERICAL RESULTS AND DISCUSSIONS 

With the theory in the above section, we proceed to calculate the Josephson current in various cases. As 
a typical case, the system temperature is taken to be zero. Besides, we take Ag = 1.0 to be the energy unit 
in this junction. 

In Figj^ we investigate the odd-FP Josephson current in Case I, and plot its spectra as a function of 
Josephson phase difference. As for the QD level and intradot Coulomb strength, we change Eq from —3.0 to 
2.0 in which U is taken to be 0.0 and 2.0, respectively. In this figure, we find that despite the change of Eq 
and U, the leading oscillation property of the Josephson current is fixed. Namely, it reaches the maximum at 
the point 9 = mr with its profile as ^ — cos9. However, the roles of QD level and Coulomb strength can 
also be clearly observed. With the departure of the QD level from energy zero point, the current amplitude 
will be suppressed gradually. Such a result is relatively apparent in Figj^a)-(b) where reflects the case of 
the zero Coulomb interaction. This can be explained as the weakness of the quantum coherence when the 
QD level is away from energy zero point. The effect of Coulomb interaction is notable in the region of Eq, 
where the QD level is occupied. It can be found that Coulomb interaction suppresses the current amplitude 
as well. This should be attributed to the destruction of the quantum coherence induced by the QD-level 
splitting (eq —7 Eq and Eq + U) in the presence of Coulomb interaction. 

According to the discussion about Case H in the second part of Sec.H, when ts gets close to tx, the s-wave 
superconductor should not be viewed as perturbation. It is natural to think that the Josephson current 
in such a case will show new properties. Next, we would like to increase the coupling strength between 
the QD and s-wave superconductor to discuss the change of Josephson effect. The numerical results are 
shown in Figj^where the tr = ts = 0.5. In this hgure, we see that in such a case, the current properties are 
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FIG. 3: Odd-FP current in Case II of [/ = 0 and U = 2.0. In (a)-(b) [7 = 0, and U = 2.0 in (c)-(d). Relevant 
parameters are chosen as It = ts = 0.5. 



FIG. 4: Josephson currents in Case III. The relevant parameters are taken to be tr = 0.1 and ts = 0.5. In the case 
of As = 1.0, the current amplitude J will be equal to be 0.04. 


completely different from those in Case I. Firstly, the current amplitude is efficiently enhanced by the increase 
of ts- Secondly, the current direction is completely reversed and its profile is deviated from Ij°^ ~ cos0. 
If the coupling between the QD and s-wave superconductor further increases, the QD will be submerged 
in the s-wave superconductor. Consider the extreme case of the weak-coupling limit (i.e.. Case III), the 
perturbation method can also be employed to evaluate the Josephson current, as displayed in the third part 
in Sec.II. It clearly shows that in such a case, the Majorana doublet couples weakly to the composite s-wave 
superconductor. Consequently, follow the relationship that = ±Jcosd and their properties are 

clearly shown [See the results in Figj^. 

In view of the current results in Case I, Case II, and Case III, one can observe that at the case of tx ^ ts 
(i.e.. Case I), the current oscillation manner is opposite to that in the other two cases. In order to clarify the 
change of Josephson current from Case I to Case III, we present the geometries of these three cases in the 
Nambu representation, as shown in Figj^ In Fig[^a)-(c), we notice that the finite coupling between the two 
MBSs of Majorana doublet, despite the direct or indirect coupling, give rise to the occurrence of anomalous 
Josephson effect. On the other hand, the coupling strength between the QD and s-wave superconductor plays 
the role in altering the inter-MBS coupling property, leading to the change of the current oscillation manner. 
In the case of tx ^ts, the MBSs couple directly to each other with a constant coupling parameter. In such 
case, the current direction is only dependent on the FP of the Majorana doublet with Ij = [j^YPJ cos 6. In 
the other case where tx ^ ts, the coupling between the QD and Majorana doublet induces the indirect 
inter-MBS coupling, as shown in Figj^a). With respect to the inter-MBS coupling in these two cases, we 
find that in the former case, the MBSs couple to each other via a nonresonant Andreev reflection process, 
whereas in the latter case, one bound state is involved in the Andreev reflection process. It is well known that 











FIG. 5: (a)-(b) Geometries of the Josephson junction in Gase I, Case II, and Case III in the Nambu representation, 
respectively. 


the quasi-particle phase will undergo a 7r-phase shift due to the presence of one bound state in the Andreev 
reflection process. Accordingly, for the case of identical FP, the current oscillations in Case I and Case III 
are opposite to each other. By the same token, we can easily see that in Case I and Case II, the oscillation 
manners of the Josephson current are opposite to each other, since an additional bound state is presented 
in the Andreev reflection process in Case II. Up to now, we have known the reason that in the considered 
junction, the current in the case of tr ^ ts is different from that in the other cases. Also, note that in 
such a structure, the role of the s-wave superconductor is to provide a channel for the coupling between the 
MBSs in the Kramers doublet and the QD is to change the channel property. For this reason, the change of 
£0 and U cannot induce any phase transition behaviors for the Josephson effect. 


V. SUMMARY 

In summary, in this work we have discussed the Josephson effects in the junction formed by a one dimen¬ 
sional Dlll-class TS and a s-wave superconductor, by embedding a QD in such a junction. Via considering 
three QD-superconductor coupling manners, we have presented a comprehensive analysis about the Joseph¬ 
son effect in this system. As a consequence, it has been found that the Josephson current oscillates in 27r 
period. However, the presence of Majorana doublet in the D///-class TS renders the Josephson current finite 
in the case of zero phase difference between the superconductors. With respect to the current direction, it is 
not only related to the FP of this junction but also depends on the coupling strength between the QD and s- 
wave superconductor. When the coupling between the QD and s-wave superconductor decreases to its weak 
limit, the direction of the Josephson current will have an opportunity to reverse. Via analyzing the particle 
motion in this structure, the reason for such a result has been clarified, namely, the QD-superconductor cou¬ 
pling manner can alter the property of the Andreev reflection between the MBSs in the Majorana doublet. 
It can be believed that this work be helpful for understanding the transport properties of the D///-class TS. 


Appendix 


According to the Hamiltonian in Eq.(19), the coupling strength between the Majorana doublet and the s- 
wave superconductor, i.e., Fo- = tt l^fcPPo-i where is the density of state in the s-wave superconductor. 
Via a straightforward derivation, one can get the result that F^ = defined by G^^{t) = 

is a retarded Green function of one QD coupled to a s-wave superconductor. By means of 
the nonequilibrium Green function technique, the matrix of the retarded Green function can be obtained, 
i.e., 


[g:; 


(w — £o)5et + iFoPo 
0 
0 

—i^^Fopo 


0 

{u> — £o)Sel + iFoPO 


0 

-i^^Fopo 
{ut + £o)Sh'f + iFopo 
0 


-i^Fopo 

0 

0 

"h £o)Sh^ + i^OpO 


(23) 




















9 


where Tq = 7rX;fc \Vks\^S{uj - ^k) and po = ■ ^<^{h)a = u 2 ±eo±u^u{ns) electron 

interaction within the Hubbard-I approximation and {rig) is the average electron occupation number. In 
the absence of magnetic factors, such a system is spin-degenerated, hence can be simplified to be 2 x 2 
matrix, i.e.. 


G 


r 

d 


{uj — eo)Se + iTopo —i^ToPo 

—i^Topo (w-I-eo)5'?i-I-iFopo 


(24) 


In the limit of strong QD-superconductor coupling, the influence of uj and £0 will be submerged, and 
can be approximated as —ipo/To, thus To- « i^po/TQ. With the increment of ts, Fq will increase. This 
surely leads to the decrease of Fo-. 
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